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Normal modes provide a fundamental basis for understanding crucial properties of solids, 
such as the thermal conductivity, the heat capacity and the sound propagationll-^. While 
the normal modes are excellently described by plane waves in crystals, they are far less un- 
derstood in disordered systems, due to the great difficulties in characterizing the heteroge- 
neous vibrational behaviors. Using charged colloids with long-range repulsion, we success- 
fully make different disordered systems without any contact friction, whose normal modes 
can be visualized at single-particle level. In these systems, we directly tackle the long-time 
outstanding puzzle in condensed matter physics: the microscopic origin of the low-frequency 
modes in disordered systems. For the first time, we experimentally clarify that the low- 
frequency modes are caused by the collective resonance of relatively disordered particles (or 
soft structures) coupled with long-wavelength transverse excitations, settling this puzzle at 
single-particle level. Next to these low-frequency modes in the density of states, we also ob- 
serve a plateau due to isostaticity, verifying the fundamental prediction of jamming model 
1 4-^. Moreover, we reveal the intrinsic correlation between the low-frequency modes and 
the real dynamics, which may lead to a universal mechanism for aging, melting and yielding. 

One important unsolved problem in condensed matter physics is to understand the normal 
modes of disordered systems. Unlike crystals, whose particles have identical surroundings and 
vibrate as plane waves, the particles of disordered systems have different local environment and 
oscillate distinctively. To fully understand their vibrations, therefore, we must correlate individual 
particles' motions with their local environment, most ideally at single-particle level. Because of 
experimental difficulties, however, such studies are rare, leaving important questions unsolved. In 
particular, disordered systems have much more low-frequency modes than crystals, which often 
form a broad boson peak in the D(uj)/uj'^^^ spectrum. These low-frequency modes are generally 
believed to cause the anomalous properties in heat capacity, thermal conductivity and acoustic be- 
haviors. The origin and properties of these modes are still open to debate. A number of models 
have been proposed. In the energy landscape model, the excess low-frequency modes are ex- 
plained as the signature from a minima-dominated phase to a saddle-point-dominated phaselvl]. In 
the jamming model, the jamming transition at zero temperature (point J) is found responsible: it 



produces ap 
cies [|4|-|6|,l8|-. 



ateau in D(uj) as the result of isostaticity, and creates extra modes at low frequen- 



1311 . Detailed numerical analysis of the low-frequency modes suggests their possi- 



ble origin as transverse phonons associated with defective soft structures fll4ll . Experimentally, 
Raman, X-ray, and neutron scatterings are commonly used to measure the normal mode spec- 
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trum. It is shown that the low-frequency modes have collective character Hi 511 and non-acoustic 
properties ni 611 ■ Most recently, the experiments in colloidal systems begin to examine the normal 
modes at single-particle level 11 17l - l20ll . Despite extensive studies, however, the microscopic origin 
of the low-frequency modes has never been experimentally clarified, which is the main focus of 
this study. 

We suspend PMMA particles of diameter a = 2.0/im (polydispersity < 3%) in density and 
refractive-index matched solvent, between two glass substrates (Methods). In the weakly polar 
solvent, theparticles and the walls carry the same type of charges and interact with long-range 
repulsion 11211 l22ll . Since there is no particle-particle and particle- wall contact, we completely 
eliminate the contact frictions in our system. Using confocal microscopy, we track every particle's 
motion with excellent spatial(15nm) and temporal (30s~^) resolutions. We study perfect 2D sam- 
ples by one single particle layer between two glass substrates. To quantitatively describe individual 
particle's local environment, we directly measure each particle's local orientational order param- 
eter, "^Qi, determined by the nearest neighbors' arrangement: \I'6j = 1 means perfect hexagonal 
arrangement around particle i and ^'e? = indicates totally random arrangement(see Methods for 
exact definition). We use different colors to represent the different \l/6j values in Fig. la. Clearly 
our heterogeneous samples contain both ordered and disordered particles, with ranging from 
to 1. We can also systematically vary the overall amount of disorder: by increasing the num- 
ber density, we increase the particle interaction and achieve more ordered samples, as shown by 
sample-A through sample-D in Fig. la and Fig. lb. 

For each sample, we track the motions of all particles for 2000 frames and extract the normal 



modes with covariance matrix method[|l8 



i(Methods). This method is valid for stable and 
metastable systems, in which no particle should rearrange; we therefore check every frame to 
make sure no particle rearranges for the entire process. Within the measurement time, every 
particle moves around its well-defined equilibrium position and the system is metastable; while for 
longer time, thermal motions cause rearrangements and the system relaxes. This relaxation limits 
our measurement time to 2000 frames. To compare the behaviors of particles under different local 
environment, in each sample we divide all particles into two groups: the relatively disordered 
particles with less than the system average (\l'6i < ^e-i), and the relatively ordered particles 
with \l/gj larger than the system average {^ei > ^ei). The comparison between the two groups 
will clarify the special role of the disordered particles in the low-frequency modes. Fig.lc shows 
the mean square displacement(MSD) for the relatively disordered particles (labeled as "dis") and 
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the relatively ordered particles (labeled as "ord") in samples A and C. The disordered curves are 
higher than the corresponding ordered ones, indicating the qualitatively larger movements of the 
disordered particles. In addition, all curves are below the dashed line of free diffusion, as the result 
of caging effect. 

From the particle tracking, we extract the normal modes at single-particle level, over three 
orders of magnitude in frequency. We plot the density of states, D{uj), in Fig.2a uM . For each 
curve, we can clearly identify two characteristic frequencies, cj^ and u^, which divide D{uj) into 
three regimes: a rapidly rising low-frequency regime {u < w^), a medium-frequency plateau 
(w^ < oj < u^), and a high-frequency regime {u > cu^). For the first time, we experimentally 
measure a plateau in D{uj), which is predicted in the jamming model as the result of isostaticity 



ll4Ll5i 



llh . This plateau did not appear in previous measurements 11 1 8l - l20h . possibly due to the 



effect of contact frictions [24]. From A to D, as system disorder decreases, the number of the low- 
frequency modes reduces, the plateau width between and uj^ shrinks, and the plateau height 
lowers. The diminishing of the plateau should come from the reduced amount of isostatic particles. 
Fig. 2b plots the D{uj)/lu^^^ spectrum and shows the boson peak in samples C, D but not in A, B. 
This is again consistent with the jamming model: the boson peak shifts to lower frequencies in 



less compressed systems such as A, B[|4|,|5L 



To illustrate the microscopic origin of the low-frequency modes, we analyze the normal modes 
at single-particle level. Fig. 3a shows three typical modes in the three different regimes of sample 
A (see Supplementary Information for other samples). The low-frequency mode is quas i-localized 
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,1130; while the 



with large-scale collective correlations, consistent with the jamming model |14|, la 
other two modes look rather random. Further inspection on the low-frequency mode reveals that 
its large-amplitude vibrations mostly occur at the disordered particle sites. To quantitatively verify 
this, we compute the average squared amplitude of the relatively disordered particles, A^iUs, and 
the relatively ordered particle, A^ord- Their ratio, A^^is/ A^oj-d, quantifies the relative vibrational 
strength between the two groups. We plot this ratio at all frequencies in Fig. 3c. In all the samples, 
A^dis/ A^ord is greater than one at small oj but quickly drops to unity as u increases. This quantita- 
tively proves that the disordered particles vibrate stronger than the ordered ones at low frequencies. 
Moreover, the transition frequency, where A^^is/ A^ord approaches unity, can separately define an 
which coincides exactly with the one found in D{ijj) (Fig. 2a). These data strongly suggest 
that all the low-frequency modes in the first regime {u < w^) are caused by the large-amplitude 
resonance of disordered particles, clarifying the microscopic origin of the low-frequency modes. 
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While the polarization vectors in real space reveal the resonance, their Fourier transform in q 
space brings further insights. In Fig. 3b, we decompose the same modes of Fig.3a into transverse 
and longitudinal Fourier components (see Methods), and show their q space distribution within the 
first Brillouin zone [— ^ < q^, qy < ^] (a is the typical distance between two neighboring parti- 
cles). We focus on the low-frequency mode in column one (a; << cu'^): the transverse component 
has much larger overall magnitude and dominates the longitudinal one, verifying the transverse 
nature of the low-frequency modes [|l4ll . Moreover, it mostly concentrates on small q values, 
indicating long- wavelength transverse excitations. At the medium frequency (w^ < tu < co^), 
the transverse component becomes uniformly distributed. When the high frequency is reached 
(to >> Lo^), the magnitude at small q decreases sharply, making the center more "black" than the 
edge. The longitudinal component keeps spreading throughout q space as to increases. 

To illustrate the overall evolution of the transverse and longitudinal components, we plot their 
average magnitude at all frequencies in Fig. 3d. For low and medium frequencies, the transverse 
curves dominate the longitudinal ones, consistent with Fig. 3b. More interestingly, each transverse 
curve also contains two characteristic frequencies, and co^, that coincide with the ones found in 
D(uj). The wonderful agreement of and co^ determined by various approaches is demonstrated 
in Fig.3e, which summarizes different properties of the normal modes from different aspects. 
Combining all the properties of the first regime (to < w^), we now draw the final conclusion on 
the low-frequency modes: they are caused by the collective resonance of the relatively disordered 
particles, coupled with long-wavelength transverse excitations. This conclusion is based on the 
direct measurements of every particle's local environment and vibrational properties in real and 
Fourier space, and may also apply to other disordered systems. 

Besides the microscopic origin, we illustrate another interesting property of the low-frequency 



modes: their ability to drive the real dynamics [|25l |26|]. By proj ecting the long-time (2000 frames) 
displacement field onto the normal mode basis (Methods) [|27ll . we find the displacement mainly 
falls on the low-frequency modes. Fig.4a and 4b compare the displacement field and the super- 
position of ten most contributed low-frequency modes in sample A, both vector fields look quite 
similar. Their difference is shown in Fig. 4c: the small residue proves their excellent agreement 
unambiguously. Fig.4d shows the sum of the projection probabilities, |C(^P, and quantitatively con- 
firms the major contributions from the low-frequency modes. The fact that we can wonderfully 
map the complex displacements of thousands of particles with only ten low-frequency modes is 
rather astonishing. It demonstrates that although the real dynamics is thermally driven, it is not 
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random during a long-time interval; instead it follows specific low-frequency modes. Similarly, 
we suspect that aging, melting and yielding may also be initiated by specific low-frequency modes, 
and speculate the excitation of these modes as the universal mechanism for these general phenom- 
ena. 

In this study, we experimentally illustrate the origin and properties of the low-frequency modes 
in disordered systems. By dividing all particles into different groups with respect to their local 
environment ("ifm), we explicitly demonstrate the important role of the relatively disordered parti- 
cles. We note that \l/6j is just one of the many parameters that can be used. For example, grouping 
particles according to their fluctuations around equilibrium positions (local Debye- Waller factor) 
can yield similar results (Supplementary Information). Therefore, our grouping approach can be 
generalized and applied to manyother disordered systems, thus providing a general method to 



analyze heterogeneous systems [|28 - 
Methods 

Sample Preparation. The PMMA particles are dyed with NBD and grafted with Polyhydroxys- 
tearic acid (PHSA) polymers. They are suspended in the mixture of lododecane and lodododecane 
(volume ratio 1:4), forming density and refractive-index matched systems. The PHSA polymers 
are slightly charged in the solvent, yielding the long-range repulsive potential between particles. 
We make 2D samples by trapping one layer of particles between two cover glasses. The glass 
surfaces are also grafted with PHSA polymers, and repel the particles from large distance. 
Local Orientational Order Parameter. The local orientational order parameter for particle i is 
defined as: '^Qi = ^ \ S^'^^e-'^^™' |. Here Ui is the number of nearest neighbors of particle i, and 
9mi is the angle between r^ — and the x axis, = 1 means the perfect hexagonal arrangement 
of six nearest neighbors and "^qi = means totally random arrangement. 

Covariance Matrix Method. From particle tracking, we can construct the covariance matrix and 
solve for the normal modes. The elements of the covariance matrix are defined as: 

= ((u.(t) - (u,,(t)))(u,(t) - (u,(t)))) (1) 

where u(t) is the particle displacement at time t, and i,j = 1, 2N run over the N particles and 
their x and y coordinates. () means time average. Diagonalizing the matrix yields eigenvectors and 
eigenvalues. Since the covariance matrix and the dynamic matrix have the relationship: Dij = 
^^{C~^)ij, the covariance matrix eigenvectors give the polarization vectors of the normal modes, 
and the covariance matrix eigenvalues, A = ^^^/''" , yield the angular frequencies uj. 



Transverse and Longitudinal Components. The transverse and longitudinal components of the 
normal modes in Fourier space are defined as: 

/u>,tr(q) = |Sj(q X ei,a;)exp(iq • r^)]^ - transverse (2) 

foj,M = |Si(q • ei,^)exp(iq • r^jl^ - longitudinal (3) 

Here u labels the mode, and are the polarization vector and the equilibrium position vector 
of particle i. The summation is over all particles. 

Projection of the Real Dynamics onto the Normal Modes. We first obtain the displacement 
field during the time window of r=2000 frames: 

\Sr{T) >= \6r{to + r) > -|5r(to) > (4) 

We then project it onto the normal mode \u! >, with the pref actor determined by: 

^/< dr{T)\6r{T) > 

The prefactors satisfy Sjct^P = 1 for a complete set of modes. We find that for some specific low- 
frequency modes, the |ca;p values are particularly large. These modes can therefore be regarded 
as the driving modes of the real dynamics. 
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FIG. 1 : Structures and motions of four samples with different amount of disorder, a, the structures by 
the particle equiUbrium positions. The total dimensions are 145 x lAb^m?. Different colors correspond to 
different ^I'gi values, b, The pair correlation function, g{r), of the four samples. The distance r is re-scaled 
by the position of the first-peak, a (1.92a, 1.89a, 1.81a and 1.73a from A to D). The peaks become more 
pronounced from A to D, indicating more order in the sample, c, MSD curves for the relatively disordered 
particles (solid symbols) and the relatively ordered particles (open symbols) in A and C. Higher MSD curves 
of the disordered particles demonstrate their larger motions. All curves are below the dashed line of free 
diffusion, as the result of caging effect. 
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FIG. 2: The distribution of tlie normal modes, a, the density of states, D{u}), of the four samples. Each 
curve can be divided into three regimes by two characteristic frequencies, and co^, as illustrated by the 
schematics. From A to D, as system disorder decreases, the number of the low-frequency modes reduces, 
the plateau width between uj"^ and oj^ shrinks, and the plateau height lowers, b, the D{oj)/oj'^^^ spectrum. 
The Boson peak appears in C, D but shifts to lower frequencies in A, B. In both a and b, the highest 
frequency plotted is determined by our spatial resolution, 15nm. 
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FIG. 3: Vibration distribution in real and Fourier space, a. The polarization vectors of three 
typical modes in sample A. The three typical modes illustrate the properties of the three regimes. 
Clearly the low-frequency mode (a; << u;^) has large-scale correlations, while the other two look 
rather random, b, The transverse and longitudinal components of the same modes in Fourier space. 
Both components are shown in the first Brillouin zone, [— ^ < qx,qy < f ]• As w increases, the 
transverse component at small q varies from very large to very small values; while the longitu- 
dinal component keeps spreading throughout q space, c, the ratio of the mean squared vibration 
amplitude between the disordered particles and ordered particles, A^dis/ A'^ord- For small u, the 
ratio is larger than one, indicating stronger vibrations of the disordered particles. The transition 
frequency where it approaches unity defines separately, d, the average magnitude of trans- 
verse and longitudinal components at all frequencies. The average is taken over the small q area 
at the center, [— ^ < qx,qy < The transverse curves (solid) dominate the longitudinal ones 
(dashed) mostly. We can again divide each transverse curve into three regimes, with two character- 
istic frequencies, and uK e, the comparison of u;^ and determined by different approaches. 
They excellently collapse onto the straight line of y — x, confirming the robustness of the three 
regimes determined differently. 
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FIG. 4: The comparison between the real dynamics and the superposition of low-frequency modes, a, 

the displacement field during 2000 frames in sample A (enlarged 2.5 times for better illustration), b, the 
superposition of ten most contributed low-frequency modes, c, the residue from a-b. The small residue 
proves excellent agreement, d, the sum of the projection probabilities, jc^jp, for all the modes. Apparently 
the low-frequency modes make major contributions. 



14 



